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Summary 

The results obtained by Kroupa, Petr & McCaughrean (1999) for specific models of young 
compact binary-rich clusters are generalised using dynamical scaling relations, to infer the 
candidate set of possible birth models leading to the Orion Nebula Cluster (ONC), of which the 
Trapezium Cluster is the core. It is found that candidate sets of solutions exist which allow the 
ONC to be in virial equilibrium, expanding or contracting. The range of possible solutions is 
quite narrow. 
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These results will serve as guidelines for future, CPU-intensive calculations of the stellar- 
dynamical and astrophysical evolution of the entire ONC. These, in turn, will be essential to 
quantify observables that will ultimately discriminate between models, thus allowing us to 
understand if the ONC is in the process of assembling a rich Galactic cluster, and, if this is the 
ON ■ case, how it occurs. 
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Subject headings: stars: binaries: general - stars: formation - open clusters and associations: 
individual (Orion Nebula Cluster, M42) 



1. INTRODUCTION 

The formation of bound star clusters is a major unsolved problem in observational and theoretical 
astrophysics. Significant observational (e.g. Lada & Lada 1991; Megeath et al. 1996; Lada, Alves & Lada 
1996; Lada, Falgarone & Evans 1997) and theoretical (e.g. Lada, Margulis & Dearborn 1984; Elmegreen & 
Efremov 1997; Klessen, Burkert & Bate 1998) efforts are being directed at solving this problem. There are 
serious technical difficulties. For example, on the observational side, the obscuration of a part of the very 
young stellar and proto-stellar population and significant uncertainties in pre-main sequence evolution (e.g. 
Wuchterl & Tscharnuter 1999) hamper the understanding of the spatial and age distribution within compact 
star-forming regions. On the theoretical side, the presently next to impossible handling of three-dimensional 
self-gravitating magneto-hydrodynamics with feedback from forming stars limits detailed understanding of 
the onset of star-formation in general, and cluster formation in particular. 

This paper forms part of a theoretical approach based on purely stellar-dynamical arguments. These are 
limited to the phase after the proto-stars and ambient gas decouple dynamically, but allow insights into the 
evolution of a stellar assemblage after the massive stars rapidly remove the remaining gas. Understanding 
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this phase of the evolution can provide a lead towards the initial conditions that give rise to bound Galactic 
clusters. 

Stellar-dynamical calculations of binary-rich young compact clusters with Aarseth's Nbody5 code 
were presented by Kroupa, Petr & McCaughrean (1999, hereinafter KPM), to study the evolution of the 
binary population and of the cluster bulk properties. The results were applied to a comparison with key 
observables available for the Trapezium Cluster, and a solution was found in which the Trapezium Cluster, 
with a stellar mass of 700 Mq, is expanding rapidly due to the expulsion of about 1000 M Q in gas roughly 
50 000 yr ago. However, while that work quantifies the evolution of the binary population and of the cluster 
within a few Myr, the question remains open as to how unique the expanding solution of the Trapezium 
Cluster is. 

In this paper the solution space of all possible models of the ONC is studied by invoking stellar- 
dynamical scaling relations that capture the gross behaviour of star clusters, and by building on the results 
obtained by KPM. Section || summarises the key observables used to constrain the allowed range of models, 
and in Section ^ the KPM results are generalised and the candidate solution space is identified. The 
conclusions follow in Section 0. 



2. OBSERVATIONAL CONSTRAINTS 

A detailed description of the observational constraints can be found in KPM. Here a brief summary is 
given. The key observables used to constrain the stellar-dynamical models by KPM are the central stellar 
number density, the pre-main-sequence age of the stellar population, the velocity dispersion, and the binary 
proportion. 

Based on a high-spatial-resolution direct-imaging near-infrared study, McCaughrean & Stauffer (1994) 
calculate that 29 systems are within the central spherical volume with a radius of R = 0.053 pc. This 
corresponds to a central number density pc° bs ' ) = 4.65 x 10 4 stars/pc 3 (strictly speaking systems/pc 3 , but 
the maximal factor of 2 error is not critical in this analysis). In a near-infrared study covering a field of 
~ 0.65 pc x 0.65 pc centred on the Trapezium Cluster, McCaughrean et al. (1996) count 700 systems. 

From the analysis of the HR-Diagram, based on optical photometry and spectroscopy, Prosser et al. 
(1994), Hillenbrand (1997) and Palla & Stahler (1999) derive a mean age for the whole Orion Nebula 
Cluster of less than about 1 Myr. 

The one-dimensional velocity dispersion within 0.41 pc of the centre of the Trapezium Cluster is 
derived by Jones & Walker (1988), from their relative proper-motion survey using photographic plates, 
to be cr( obs ) = 2.54 ± 0.27 km/s. Because the plate-reduction algorithms eliminate any signature due to 
rotation and/or expansion or contraction, it is unknown from that work and others (see KPM) if the 
Trapezium Cluster is expanding or contracting. Frink, Kroupa & Roser (1999) measure the bulk motions in 
the ONC using new absolute proper motion data, and find some expanding as well as contracting motion. 
Their result is consistent with the ONC currently expanding in an approximately cylindrical molecular 
cloud potential (Kroupa & Frink 1999). Following KPM the velocity dispersion obtained from the more 
precise relative proper motions is used in the present analysis. 

High-resolution imaging using a variety of techniques has shown that the binary proportion within 
about 0.3 pc of the centre of the Trapezium Cluster is significantly below that seen in Taurus-Auriga, but 
consistent with the Galactic field value (Prosser et al. 1994; Petr 1998; Petr et al. 1998). This may be the 
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result of disruption through encounters in the cluster (KPM). 

The Trapezium Cluster contains virtually no gas (Wilson et al. 1997), and a large part of the ONC 
population is optically visible (Hillenbrand & Hartmann 1998), so that purely stellar-dynamical calculations 
can be applied. 

There is little evidence for sub-structure in the ONC (Bate, Clarke & McCaughrean 1998). 

3. THE KPM STELLAR-DYNAMICAL MODELS 

Details are found in KPM. The KPM models contain A = 1600 stars (cluster mass Ma = 700 Mq), with a 
Galactic-field IMF with stellar masses in the range 0.08 M Q to 30 Mq. The stars are paired at random to 
form a primordial binary proportion of unity, / to t = 1 with a period distribution as seen in Taurus-Auriga, 
and /tot = 0.6 with a period distribution as in the Galactic field. Here, the total (counting all separations) 
binary proportion is /tot = A D in/(Abi n + A s ; ng ), where Abj n and A s ; ng are the number of binary and 
single-star systems, respectively. Stellar evolution is not treated, and the dynamical evolution is followed 
for 5 Myr only. 

Six different models are constructed to cover a range of possible cases, each having initially a Plummer 
density profile with half-mass radius Rq.5- 

• Virial equilibrium models with primordial /tot = 1 (model Al) and /tot = 0.6 (model A2). These have 
initially Rq, 5 = 0.1 pc, i cross = 0.06 Myr (eqn. ^) and t ro i ax = 0.62 Myr (eqn. |). 

• Expanding clusters assuming a star formation efficiency e = Ma/ (Ma + M gas ) — 0.42 and 
instantaneous loss of a gas mass M gas = 967 Mq at the start of the computation, with binary 
proportions /tot = 1 (model Bl) and 0.6 (model B2). These models have i?o.5 = 0.1 pc and an initial 
velocity dispersion that is e -1 / 2 larger than the virial equilibrium value without the gas. 

• Clusters that collapse with an initial velocity dispersion that is a factor 0.14 times smaller than the 
virial equilibrium value, with initially i?o.5 = 0.4 pc (model CI) and 0.8 pc (model C2). These have 
/tot = 1. 

Three renditions of each model yield averages of the relevant observables. 

The choice for A is based on the observed number of systems in the Trapezium Cluster. The initial 
central density (p c ,in = 10 592 stars/pc 3 ) of the virial equilibrium and expanding models is larger than the 
presently observed value (/9 c ° bs ' ) = 10 4 67 stars/pc 3 ), because a reduction of the central density with time 
is anticipated without knowing the detailed evolution in advance, since no such compact and binary-rich 
clusters had ever been evolved until then. For the collapsing models, the initial possible configuration is 
unknown, since A-body collapse calculations of an initially centrally concentrated density profile and with 
a large proportion of realistic binary systems had never been done before. The values i?o.5 = 0-4 and 0.8 pc 
are thus first tries. 

The mean stellar mass is 0.44 Mq. It is used throughout this paper: Ma = 0.44 x N Mq. 
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4. GENERALISATION 

Before constraining the solution space using dynamical arguments, a pre-selection can be made immediately, 
by noting that the lower limit on N is 700 stars, since this many systems have been counted in the 
Trapezium Cluster (Section |||). 

It is known that the binary fraction is high (about 50-60 per cent), so N is likely to be closer to 1200. 
However, since some loss of stars is expected through dynamical ejection, and in general location at larger 
radii than the survey limit, N = 1600 is deemed more appropriate to be a lower limit. There is thus little 
room for an initial N below about 1600. 

Assuming the Trapezium Cluster is merely the core of the ONC, an upper limit N — 5600 — 10000 is 
obtained for the whole ONC. The lower value is valid if the overall binary proportion is /tot = 0.6, since 
about 3500 systems have been counted within approximately 2.5 pc of the cluster's centre (Hillenbrand & 
Hartmann 1998). The upper limit is estimated by assuming / tot = 1 as a probably more realistic value for 
most of the ONC, implying that 7000 stars will have been counted, and that 30 per cent of all systems have 
been lost through expansion beyond the 2.5 pc radius, or have not been seen because they are too deeply 
embedded or too faint. 



4.1. Clusters in virial equilibrium 

Without doing additional time-consuming numerical experiments, the evolution of model clusters in virial 
equilibrium with different N and i?o.5 can be assessed through simple scaling arguments. 

Binary depletion occurs approximately on a crossing time-scale, Across (fig- 3 in KPM), whereas the 
velocity dispersion and central density change on a time-scale comparable to the relaxation time, i re iax (fig- 2 
in KPM). That complete destruction of long-period or soft binaries occurs within a few tens of initial £ C ross 
is demonstrated in Fig. [l|, in which the evolution of /tot for two very different clusters is compared. Despite 
very different initial crossing times and different N, /tot decays on a comparable time-scale ~ few x t c 
That the decay is not exactly in-phase does not affect the conclusions of this paper. 
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The three-dimensional velocity dispersion, 

'GM, 
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where G = 4.49 x 10~ 3 pc 3 /(Mo Myr 2 ) is the gravitational constant, and s is a structure constant of 
order unity defined by the details of the density distribution. The one-dimensional velocity dispersion 
is o id = a/y/3, assuming the velocity distribution is isotropic. For M c \ = 700 M Q and R0.5 = 0.1 pc, 
0id/s = 2.29 pc/Myr (1 pc/Myr= 0.982 km/s, but throughout this paper no distinction is made between 
these two units), and from fig. 2 in KPM, <rm(t = 0) = 1.95 km/s, so that s — 0.85 throughout this paper. 



The typical crossing time, 
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whereas the relaxation time in Myr (from Binney & Tremaine 1987), 
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Thus, it is immediately apparent that increasing N (i.e. M c \) will, for the same R0.5, shorten the 
time-scale for binary depletion while simultaneously slowing down the time-scale for the evolution of the 
central density and velocity dispersion. This implies that by increasing N, the variation of the central 
density becomes increasingly negligible over the age of the cluster. This is already evident in fig. 2 in KPM, 
where the central density decays by a factor of 2.5 only during the first Myr. The central density for a 
Plummer density distribution is 

Pc = 5: (4) 

4ir (0.77i? . 5 ) 3 

where R p \ = 0.77 i?o. 5 is the Plummer radius. 

For illustrative purposes, two sets of initial models are constructed and compared in Fig. |^. One 
set of models has log 10 (p c ) = 5.92 (stars/pc 3 ), which is the same as in models A. The other set has 
log 10 (p c ) = 4.67 (stars/pc 3 ), which is the same as the observed central density in the Trapezium Cluster. 
The first set is clearly not a viable solution set since the initial density remains too high (t re iax is too long). 

However, there is a set of virial-equilibrium candidate solutions if 10000 A < 16000 and 
0.48 ^ R0.5 ;$ 0.61, for the constraint p c = constant. These models have a velocity dispersion consistent with 
the observed value, and t CTOSS = 2.5 x 10 5 yr, which is sufficiently short for significant binary depletion to 
have occurred in the ONC until now (Fig. [[]). Fully self-consistent A-body calculations are needed to verify 
if a Taurus-Auriga-like binary population can be reduced to the observed level in the Trapezium Cluster. 
The calculations of KPM suggest that this is likely to be the case. 

Relaxing the constant-central density constraint, the full set of candidate solutions is obtained (Fig. ||). 
As already apparent in Fig. ||, the set of presumed solutions is very narrow, with A w 10 4 , M c \ » 4.4 x 10 3 Mq 
and i?o.5 ~ 0.45 pc. It has been known for a long time that the velocity dispersion in the Trapezium Cluster 
is super-virial, if only the number of detected stars is used to estimate the cluster mass (e.g. Hillenbrand & 
Hartmann 1998). However, the candidate solutions found here demonstrate that, by taking into account the 
binary companions, the natural loss of stars from the cluster through its dynamical evolution, and allowing 
for a small fraction of stars not detected yet, a model of the initial ONC is found that is remarkably close 
to virial equilibrium. 



4.2. Expanding clusters 

Model B2 in KPM shows that if the Trapezium Cluster is expanding then it is being observed now at a time 
when the velocity dispersion and central density have not decayed so much as to make the system appear 
like an association. Expansion must have begun a few tens of thousands of years ago, and the initial / to t 
must have been similar to that in the Galactic field, and thus significantly below that in Taurus-Auriga. 
Model B2 assumes a high star-formation efficiency of e = 0.42, which nevertheless leads to an essentially 
rapid free expansion of the entire system (Kroupa & Frink 1999). Indeed, Lada, Margulis & Dearborn 
(1984) demonstrate that, if e < 0.5 and gas is removed instantaneously, then an unbound, expanding 
association results. The presence of a few O stars in the Trapezium Cluster suggests that the removal of 
gas occurred on a time-scale comparable to or shorter than a crossing time. 

The set of candidate solutions to an expanding Trapezium Cluster, assuming instantaneous gas removal 
and virial equilibrium of the stellar and gaseous system before gas expulsion, can be constrained by noting 
that the pre-expansion central density and velocity dispersion must have been larger than the presently 
observed values, p c (t = 0) > pi ohs ^ and axo(t — 0) > a^ ohs \ where t — refers to the gas-expulsion time. 
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The central density cannot be made arbitrarily large, since binary systems with semi-major axes as large as 
a = 440 AU exist in abundance (section 2.3 in KPM). Thus, if each binary system needs a volume with a 
radius 7 a to form, and the forming systems touch, then the maximum number density would be 

max _ 3X2 

so that p™ ax = 10 7 - 69 - 10 4 - 69 stars pc~ 3 for 7 = 1 - 10. 

In a freely expanding flow that is unhindered by self-gravity, the velocity dispersion measured within a 
projected volume with a constant radius 2 R decreases with time according to 

* 1D(t) = f + a/L(0)' (6) 



where <7id(0) = owyj^fl is the velocity dispersion before gas expulsion, and ctid is, as in Section 4.1, the 
velocity dispersion of the stellar cluster if it were in virial equilibrium and without gas (eqn. Q). This 
equation ensures that the initial velocity dispersion, ctid(0), is arrived at, and essentially measures the 
typical velocity, R/t, of the stars remaining in the measurement area. The term R/am(Q) is the time-lag 
until stars with typical velocity 0id(O) leave the area. The velocity-dispersion model-data corrected for the 
radial bulk flow plotted in the bottom panel of fig. 6 in KPM fit this relation very well with R — 0.21 pc. 
That the flow in the KPM models B is free on scales larger than a few tenths of a pc is also demonstrated 
in Kroupa & Frink (1999). The time at which oi D (i exp ) = er( obs ) = 2.54 ± 0.27 km/s is 

In a free radial expansion flow, the central density in a Plummer density distribution should evolve as 

because the Plummer radius of the system increases with time as R v \{t) — 0.77 Rq. 5 + a(0)t. However, the 
model data plotted in the upper panel of fig. 6 in KPM show that the real central density, evaluated within 
a radius R — 0.053 pc, decreases much less rapidly. This is due to stellar-dynamical interactions in the 
binary-rich system leading to a part of the population decoupling from the flow and forming a small bound 
cluster, as discussed in section 5.2.1 in KPM. Cooling of the flow through some binary destruction appears 
to be instrumental in forming this decoupled core, as can be inferred by comparing the central densities 
of models Bl and B2 at t > 0.2 Myr in fig. 6 of KPM. This mode of cluster formation will be addressed 
in more detail in a future contribution. For the present, it has to be accepted that the evolution of p c {t) 
cannot be expressed analytically. 

A set of possible initial models is summarised in Fig. |[ Surmised solutions with N = 1600 and a 
star-formation efficiency of 42 per cent have R0.5 < 0.12 pc, ensuring that the pre-expansion velocity 
dispersion is larger than the observed value by at least a one-sigma error margin. Such models lead to the 
velocity dispersion being consistent with the observed value within 10 5 yr. For e = 0.1, the candidate 
solutions have R0.5 < 0.27 pc, ensuring the pre-expansion central density is larger than the currently 
observed value. In this case the pre-expansion velocity dispersion is very large, requiring a longer time 
to pass until the velocity dispersion has decayed to a value consistent with that observed. Nevertheless, 
i oxp < 10 5 yr. Candidate solutions with N = 10 4 have i? 05 < 0.47 pc and t cxp < 10 5 yr for e < 0.42. 
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The full set of candidate solutions is shown in Fig. |5|. For all candidate solutions, i cxp < 1.2 x 10 5 yr. 
The KPM calculations show that /tot does not evolve significantly, so that /tot ~ 0.6 is required at the time 
when expansion ensues. 

4.3. Collapsing clusters 

An initially homogeneous sphere with radius i?i n that starts from rest collapses to a singularity under 
its own gravity within a free-fall time, i co i ~ t CTOSS /^/2 (e.g. Binney & Tremaine 1987). If the initial 
matter distribution is not smooth, for example when only a finite number of stars is present, local density 
differences will lead to the growth of the tangential velocity dispersion which limits the radius of the collapse 
to the minimum value at the bounce (Aarseth, Lin & Papaloizou 1988) 

R h « i? in /iV 1/3 . (9) 

The post-collapse system exhibits a core-halo structure, and an initial asphcricity will be retained to some 
degree by the halo, whereas the core will be approximately spherical (Boily, Clarke & Murray 1999). This 
may, in principle, be a model of the ONC, which has an approximately spherical core (the Trapezium 
Cluster) and an elliptical halo (Hillenbrand & Hartmann 1998), provided the ONC is old enough for the 
collapse to have progressed to the morphology evident today. 

The presence of a large number of primordial binary systems complicates the situation dramatically 
because they have relatively large interaction cross sections, and they constitute an energy sink as well as a 
highly significant energy source. A cold collapse can be envisioned being arrested when the densities become 
large enough that the binary systems begin interacting. At this stage cooling through binary disruption 
will reduce the relative velocities between the systems, slowing down the collapse. Additionally, heating 
through hardening binaries will also oppose the collapse. Vesperini & Chernoff (1996) address such issues 
by considering the collapse of an initially homogeneous sphere with a binary proportion /f ot = 0.05. They 
find that the primordial binaries do not alter the collapse and violent relaxation. However, their results are 
not applicable to the present work, since the binary proportion used here is significantly higher, and because 
the results of Vesperini & Chernoff apply to the case where binary-binary collisions are insignificant. 
Binary-binary collisions, however, dominate the interactions in the realistic case considered here. Details of 
the processes await further research. 

During a cold collapse starting from an initially spherical and homogeneous configuration with initial 
(central) density p c : m , the central density increases to a value at the bounce p c ^ N x p c ,m, given the 
above scaling for the linear dimension of the system. This is not fulfilled in models CI and C2, since this 
scaling and eqn. || for p c ,in, would imply p^\ — 10 7 3 [stars/pc 3 ] and p®\ = 10 6 5 [stars/pc 3 ], whereas from 
fig. 10 in KPM values of 10 5 9 and 10 5 8 stars/pc 3 , respectively, result in the collapse computations. The 
central density at bounce is significantly smaller because (i) the models have a small but finite initial 
velocity dispersion, (ii) the initial density distribution is centrally concentrated leading to a spread of arrival 
times near the centre which increases i?b, and (iii) binary systems oppose the collapse. 

Similarly, the one-dimensional velocity dispersion of the system at bounce cannot be calculated from a 
simple scaling law, because at this stage it is unknown how the binding energy of the binary-rich initially 
centrally concentrated cluster evolves. However, the data in the lower panel of fig. 10 in KPM suggest 
that, if £r v i ri j n is the one-dimensional velocity dispersion a pre-collapse cluster would have were it in virial 
equilibrium (eqn. |l|), then 

(7b = 1.71 X OViran- (10) 
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For models CI and C2, the one-dimensional velocity dispersion er v ir,in = 0.97 km/s and 0.69 km/s (eqn. [j]), 
respectively, so that a^ 1 « 1.62 km/s and o^ 1 « 1.20 km/s. Useful constraints can be placed on the 
candidate solution-set of pre-collapse models of the ONC if eqn. [l(] is adopted for a comparison with cr^ obs ' . 

The one-dimensional velocity dispersion at bounce must be larger than or equal to the observed value 
<7b ^ cr( obs ). Also, the initial central density must be smaller than or equal to the presently observed value, 
Pc,in < Pc° bs ^- In addition, the collapse time must be short enough to be consistent with the pre- main 
sequence age of the Trapezium Cluster, t co \ ^ 1 Myr, since the cluster could not have been collapsing 
stellar-dynamically for a time longer than the age of the stars. 

The collapse time t co \ oc i C ross, where t cross is the crossing time of the pre-collapse cluster (eqn. ||). To 
estimate the constant of proportionality, the data presented in fig. 10 of KPM need to be resorted to again. 
The result is 

t co \ = 0.97 x t cross . (11) 

Sequences of initial models are constructed for A — 1600 and 10 4 , and p c ,m, tool and Ob are calculated 
for each as a function of i?o.5- The results are presented in Fig. ^. There is no candidate set of solutions for 
N = 1600, because the velocity dispersion at the bounce never reaches the observed value for those cases 
when the central pre-collapse density is smaller than the observed central density in the Trapezium Cluster. 
However, there is a candidate solution set for models with N — 10 4 . Such models have 0.5 < i?o.5 < 1-2 pc 
and should satisfy the observational constraints at some time ^ 1 Myr during collapse. Detailed TV-body 
calculations are required to verify if this is indeed a viable solution set, and if the central density and 
velocity dispersion pass through the observational constraints at about the same time. 

The full set of candidate solutions is shown in Fig. 0. The results of KPM illustrate that /tot decreases 
significantly during 1 — 2 x t co \, but that the binary proportion after virialisation depends on the initial 
Rq.5- Thus, it is likely that initial / tot < 1 (i.e. a smaller binary proportion than in Taurus- Auriga) will be 
necessary for presumed solutions with Rq^ > 0.5 pc, to account for the observed binary proportion in the 
Trapezium Cluster. Detailed A-body calculations are needed to address this issue. 

4.4. Caveats 

The models investigated here assume the clusters initially have spherical and smooth Plummer density 
profiles. Deviations from this will not change the results significantly for the same p c , M c \ and i?o.5- A 
complication that will have to be incorporated in future studies of the evolution of the ONC is the tidal 
field exerted across it from the nearby molecular cloud. Kroupa & Frink (1999) show that the overall shape 
and velocity field of an expanding young cluster can be affected noticeably in such a case. 

5. CONCLUSIONS 

Using stellar-dynamical scaling laws together with results from A-body calculations, three general classes 
of possible stellar-dynamical solutions to initial configurations of the ONC (and of the Trapezium Cluster) 
are investigated. By doing so, candidate initial models have been found for all three classes. Thus, the 
ONC may presently be (i) in or close to virial equilibrium, (ii) expanding, or (iii) collapsing or in violent 
relaxation following a cold collapse. Class (i) only allows a candidate solution in virial equilibrium with 
Ro.5 ~ 0.45 pc and A w 10 4 . Such a large A is realistic only if the binary proportion in the ONC (rather 
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than the Trapezium Cluster) is close to unity, and if about 30 per cent of stars are not detected and/or 
have already been ejected to radii of a few pc or larger. 

The first and third classes will undoubtedly lead to a bound Galactic cluster, while the second is likely 
to, on the basis of the formation of a small bound cluster in the expanding Trapezium-Cluster models 
of KPM. The present results thus show that the ONC is most probably forming a Galactic cluster. It is 
unclear though by which of the three paths above this is occurring. 

However, only candidate initial models have been presented here. Follow-up, CPU-intensive iV-body 
calculations will have to be performed to verify which of these are consistent with the observational 
constraints. The constraints include the velocity dispersion, mass segregation, and the binary-proportion 
and density profiles in the ONC. 
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Fig. 1. — Depletion of the binary proportion in two very different clusters. The time-scale is in units of the 
initial crossing time, whereas the vertical axis is scaled to the minimum and maximum /tot- The long-dashed 
curve (I) is the average of five iV-body renditions of a cluster with initially 200 binaries and a Plummcr density 
distribution with i? .5 = 2.5 pc (p c = 13 stars/pc 3 ), and stars with masses in the range 0.1 to 1.1 M Q (from 
Kroupa 1995). The initial crossing time is 18 Myr, the initial / to t = 1.00, and the final /tot = 0.83. The 
short-dash-dotted curve (II) is for model Al in KPM (and this paper; p c = 10 5 ' 92 stars/pc 3 ), for which 
the initial / to t = 0.76 (reduced from unity because of disruption through crowding, see KPM), the final 
/tot = 0.34, and initial t CIOSS = 0.059 Myr. For the virial-equilibrium models discussed in this paper, /tot (*) 
will lie between the two extreme cases shown here. 
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Fig. 2.— Two sets of initial model clusters in virial equilibrium indicated by the super-scripts M (long- 
dashed curves) and O (short-dash-dotted curves). M refers to models that have the same initial density 
as models A, and O refers to models with a central density equal to the observed value in the Trapezium 
Cluster. The number of stars, 77 = N, the relaxation time, r] = i ro i ax , and crossing time, r\ — t cross , and 
the one-dimensional velocity dispersion, 77 = ctid, are plotted as logio[??(i?o.5)] for each model. For example, 
models A have initially -R0.5 = 0.1 pc with N = 1600, logio(t C ross) = —1-23 (Myr), logio(irciax) = —0.21 
(Myr) and logio(criD) = 0.29 (km/s). These are identified by the large black dots. The thin vertical dotted 
lines delineate the possible set of equilibrium solutions: O-models with 0.48 < i?o.5 < 0.61 pc have the 
same central density as the Trapezium Cluster, and a velocity dispersion consistent within the one-sigma 
uncertainty range with the observed value, which is shown by the two horizontal dotted lines. 
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Fig. 3. The full set of candidate initial cluster models in virial equilibrium. The solution set is obtained by 
constructing a grid of models with different i?o.5 and N, and selecting only those that simultaneously satisfy 
2.27 < c id < 2.81 km/s, 23.6 < n c < 34.4 stars and i C ross < 0.50 Myr, where n c is the number of stars in 
the central volume with radius R — 0.053 pc. The boundaries on ctid are the one-sigma uncertainty range 
on the observed value. The limits on n c are 29 ± -s/29- Note that N > 10 4 must also be excluded because 
this many stars were certainly not present in the ONC when it was born. 



-14- 




U (e=0.42) 







-J I L. 



-J I L_ 







0.2 



0.4 
R 



_i i i_ 



-J I L. 



— 1 I L_ 



0.5 



0.6 
[pc] 



0.8 



Fig. 4. — Three sets of possible initial model clusters that expand as a result of instantaneous gas loss. 
The upper panel displays two sets of models for N — 1600, and the lower panel shows one set of models 
for N = 10 4 . For each model the initial central density, logio(p c ), is plotted as the solid line, the prc- 
expansion velocity dispersion, ctid(O), is plotted as thick short-dashed curves for star-formation efficiencies 
e = 0.42 and 0.10, and the time, logio(t e xp), when ciD^cxp) = 1-73 km/s (the three-sigma lower limit on 
the observed velocity dispersion) is shown as the dash-dotted line. The lower three-sigma value is used 
to provide an upper boundary on t cxp for each i?o.5- The one-sigma uncertainty range on the observed 
velocity dispersion is represented by the two lower horizontal dotted lines, and the upper horizontal dotted 
line is the observed central density, logio(/9c° bs ' ) ) in the Trapezium Cluster. For example, for model B2, 
logio(Pc) = 5.92 [stars/pc 3 ], cr 1D (0) = 3.01 km/s and log 10 (i oxp ) = 0.72 [10 4 yr], i.e. t cxp = 5.2 x 10 4 yr 
(large black dots). This is the solution found by KPM. 
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Fig. 5. — The full set of candidate initial models for expanding clusters. The candidate solutions satisfy 
simultaneously ctid > 2.81 km/s, n c > 34.4 stars (see Fig. ||) and p c < 10 7 ' 69 stars/pc 3 . Black dots are for 
a star-formation efficiency e = 0.42, and open squares are for e = 0.1. Note that logio(-/V) < 3.5 should be 
excluded because at least 3500 ONC stars have been seen, and N > 10 4 must be excluded because this many 
stars were certainly not present in the ONC when it was born. 
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Fig. 6. — Two sets of possible pre-collapse model clusters. Long-dashed curves are for TV = 1600 and dash- 
dotted curves are for ./V = 10 4 . For each, logio(?7) i s plotted in dependence of R0.5, with r\ = p c ,in,t C oi,o'b- 
For example, model CI has p c ; m = 1.3 x 10 4 stars/pc 3 , and from fig. 10 in KPM, t co \ — 0.45 Myr and 
(Tb = 1.62 km/s (black squares), and model C2 has p Ci ; n = 1.6 x 10 3 stars/pc 3 , and from fig. 10 in KPM 
t co \ = 1.32 Myr and (Tb = 1.20 km/s (black circles). The upper horizontal dotted line is the presently 
observed central density in the Trapezium Cluster, and the two lower horizontal lines show the one-sigma 
range of the observed velocity dispersion. The two vertical dotted lines delineate the candidate solution set 
for N — 10 4 , the left boundary being determined by p c ^ n < pi ^ and the right boundary by Cb > cr' obs ^. 
Models with 0.5 < R0.5 < 1.23 pc are thus expected to satisfy the observational constraints as they collapse 
with collapse times t co \ < 10 6 yr. There is no set of solutions for models with N = 1600, because when 
(7 b > a(° hs \ /9 Cj ; n is larger than the observed density. 
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Fig. 7. — The full set of candidate initial models of collapsing clusters. The candidate solutions satisfy 
simultaneously (Tb > 2.27 km/s, n c < 23.6 stars (see Fig. ||) and t co \ < 1 Myr. Note that logio(iV) < 3.5 
should be excluded because at least 3500 ONC stars have been seen, and N > 10 4 must be excluded because 
this many stars were certainly not present in the ONC when it was born. 



